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Abstract 


Sensors on Landsat satellites have been collecting images of the Earth’s surface for nearly 40 


years. These images have been invaluable for characterizing and detecting changes in the land 
cover and land use of the world. Although initially conceived as primarily picture generating 
sensors, even the early sensors were radiometrically calibrated and spectrally characterized prior 
to launch and incorporated some capabilities to monitor their radiometric calibration once on 
orbit. Recently, as the focus of studies has shifted to monitoring Earth surface parameters over 
significant periods of time, serious attention has been focused toward bringing the data from all 
these sensors onto a common radiometric scale over this 40- year period. This effort started with 
the most recent systems and then was extended back in time. Landsat-7 ETM+, the best- 
characterized sensor of the series prior to launch and once on orbit, and the most stable system 
to date, was chosen to serve as the reference. The Landsat-7 project was the first of the series to 
build an image assessment system into its ground system, allowing systematic characterization 
of its sensors and data. Second, the Landsat-5 TM (still operating at the time of the Landsat-7 
launch and continues to operate) calibration history was reconstructed based on its internal 
calibrator, vicarious calibrations, pseudo-invariant sites and a tie to Landsat-7 ETM+ at the time 
of the commissioning of Landsat-7. This process was performed in two iterations: the earlier 
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one relied primarily on the TM internal calibrator. When this was found to have some 
deficiencies, a revised calibration was based more on pseudo-invariant sites, though the internal 
calibrator was still used to establish the short-term variations in response due to icing build up 
on the cold focal plane. As time progressed, a capability to monitor the Landsat-5 TM was 
added to the image assessment system. The Landsat-4 TM, which operated from 1982-1992, 
was the third system to which the radiometric scale was extended. The limited and broken use 
of the Landsat-4 TM made this analysis more difficult. Eight-day separated image pairs from 
Landsat-5 combined with analysis of pseudo invariant sites established this history. The fourth 
and most challenging effort was making the Landsat-1 to -5 MSS sensors’ data internally 
radiometrically consistent. This effort was particularly complicated by the age of the MSS data, 
varying formats and processing levels in the archive, limited datasets, and limited 
documentation available. Ultimately, pseudo-invariant sites were identified in North America 
and used for this effort. Note that most of the Landsat-MSS archived data had already been 
calibrated using the MSS internal calibrators, so this processing was imbedded in the result. The 
final effort was developing an absolute scale for Landsat MSS similar to what was already 
established for the “TM” sensors. This was accomplished by using simultaneous data from 
Landsat-5 MSS and Landsat-5 TM, accounting for spectral differences between the sensors 
using EO-1 Hyperion data. The recalibrated history of the Landsat data and implications to 
users are discussed. The key result from this work is a consistently calibrated Landsat data 
archive that spans nearly 40 years with total uncertainties on the order of 10% or less for most 
sensors and bands. 
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1. Introduction 

1.1 Overview 

With the launch of Landsat-1 in 1972, a new era in repetitive global coverage of the Earth’s 
surface with moderate spatial resolution imagery commenced. Initially, much Landsat data 
analysis was by visual interpretation, or if by computer, involved training the computer to 
recognize certain targets based on within-scene statistics, i.e., classification. Although the 
Multi- Spectral Scanner (MSS) sensors on these early Landsat satellites were spectrally 
characterized and radiometrically calibrated prior to launch and did include capabilities for on- 
orbit radiometric calibration, there was little demand for consistent quantitative radiometric data 
and the emphasis in data processing was on limiting artifacts, a.k.a., destriping and producing 
visually consistent products. With the transition to the Thematic Mapper (TM) sensors in the 
early 1980s, the value of radiometrically calibrated data began to be recognized, i.e., retrieving 
Earth surface parameters by converting the sensor output to physical units (radiance) and then 
atmospherically correcting the data. These early analyses were severely hampered by the high 
cost of Landsat data products (up to $4000, at times), resulting in part from the Landsat 
commercialization effort, particularly for time series analyses that required multiple Landsat 
scenes. 

With the return of the Landsat system to US government control beginning in 1992 with the 
development Landsat-7 with its Enhanced Thematic Mapper Plus (ETM+) sensor and later with 
Landsat-4 and -5 and their data archives, the data costs eventually were on the order of 
$600/scene. This increased demand for these products. At the same time, the US government 
(jointly between NASA and USGS) enhanced the Landsat ground processing system to include 
a staffed Image Assessment System (IAS) which provided for routine characterization of the 
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geometry and radiometry of the data and calibration updates as necessary, significantly 
improving the quality and consistency of the data products. When the Landsat-5 system reverted 
to USGS control, efforts to bring these data into harmony with the recently launched and 
calibrated Landsat-7 ETM+ peaked and the IAS was expanded to include Landsat-5 TM data. 
Since that point in time there has been a steady effort to bring all the Landsat systems into a 
consistent non-proprietary processing system that could be readily updated to provide products 
over the lifetime of Landsat. Dropping all user costs for the US Landsat data products, 
beginning in 2008, has tremendously increased the demand for the data and has further driven 
this effort. The approach has been to systematically work backward from the present (Landsat- 
7 ETM+ and Landsat-5 TM) to the past. Papers have been published about the details of much 
of this effort over the past 10 years or so and a number of papers are in the process of being 
written on the last pieces of this effort. This paper will be summarizing analyses, results and 
implications for data users, and will rely heavily on referencing other work for the details of the 
process. 

1.2 History of Satellites and Sensors 

Table 1 provides a summary of the history of the Landsat series of satellites, including dates of 
operation and sensors (NASA, 2011). A power problem related to the Return Beam Vidicon 
(RBV) system stopped its use after the first two weeks on orbit for Landsat- 1 and although the 
RBV functioned on Landsat-2, it was used little. The RBV was modified for Landsat-3, 
converting it to a single band, higher spatial resolution (circa 30 meter) system. The RBV will 
not be discussed further in this paper. The Landsat-3 MSS had a thermal band, but it did not 
perform well. The TM’s and ETM+ also have a thermal band, whose calibration history is 
discussed in a companion paper (Schott et al., 201 1). 
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1.3 Basic Sensor Descriptions and Radiometric Calibration Capabilities 

Table 2 provides a summary of the MSS, TM and ETM+ characteristics for the reflective 
portion of the solar spectrum. The MSS on Landsats 4 and 5 was modified to accommodate the 
lower altitude of Landsat-4 versus the previous Landsats. Also, its solar calibrator was 
removed. 


2. Evolution of Calibration Record 

2.1 State as of 1999 

During the time of the development and up to the launch of Landsat-7 in April 1999, the data 
products available in the United States from the currently operating and historical Landsat 
sensors were in a variety of formats from several different processing systems. Processing 
systems were operated by the United States Geological Survey (USGS) and by a commercial 
operator. Some data products for US government users were produced on a different system 
from those available to the general public. Most of the processing systems were proprietary 
and publicly available information was limited. The insight into how these processing systems 
operated was rather limited, even by the USGS as their operator; and the USGS’s ability to 
update these systems with new constants or algorithms was complicated and not routinely 
performed. In general, data were processed with the same calibration constants and algorithms 
in place as when the satellites completed their missions and these constants and algorithms 
could be a function of both acquisition data and processing date. The state of the MSS 
radiometric processing systems is partially documented in Markham et al, (1985), Helder 
(1993), Markham and Barker (1986) including the radiometric algorithms being used and the 
radiometric scaling parameters for the final products. 
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2.2 Establishment of Landsat-7 ETM+ as the standard: circa 2001 


2.2.1 Landsat-7 mission improvements 

A number of factors converged during the development of the Landsat-7 sensor and mission 
that led to a significant step forward in the quantitative use the sensor data. Some of these were 
a result of improvements in technology in the 15 or so years between when the TM’s and the 
ETM+ were designed and fabricated. Others were a result of “lessons learned” from the earlier 
programs and the involvement of some core personnel on the instrument teams, science team 
and calibration/validation teams that had previous Landsat TM experience. Finally, some were a 
result of an evolution of the science data user community including the greater recognition of 
the value of continued calibration and characterization efforts for satellite sensors throughout 
their missions. These led to an improved sensor in terms of stability (both gain and offset), 
freedom from some image artifacts e.g., memory effect, scan correlated shift, (Helder et al, 
2004), absolute calibration and characterization both pre- and post-launch, and the 
comprehensiveness of the acquired data. 

Specific examples include: 

A. Flight hardware 

1. Ion- Assisted Deposition (IAD) significantly improved the stability of the spectral 
filters, both between ambient and vacuum conditions and over time in general. 

2. A/D convertors were markedly improved between F5 TM and F7 ETM+, particularly 
in the area of step size uniformity. 
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3. Contamination control increased in visibility within NASA programs. The control of 
contamination on instrument and optics likely contributed to minimizing ETM+ optic 
throughput degradation 

4. ETM+ design changes as a result of L4/L5 TM issues. Electronics changes 
effectively removed the memory effect and scan correlated shifts. The addition of an 
optical window in front the cold focal plane eliminated the “ice” buildup on the cold 
focal plane and resulting fluctuations and degradations in cold focal plane band 
responsivity. 

5. The addition of on-board calibration devices: (a) solar diffuser, that generally worked 
well, though it degraded, and (b) partial aperture solar calibrator that turned out to be not 
very useful. 

B. Science Data Processing. 

1. Non-propriety and publicly released data processing algorithms and software 
improved the understanding and ability to modify parameters and algorithms for 
radiometric and geometric processing. 

2. A staffed image assessment system (IAS) was built into the data processing system. 
The IAS is an off line tool that is used to characterize sensor performance and data 
quality and update parameters regularly and algorithms as required. 

C. Landsat Science Team and Land Cover Project Science Office 

1. Though short lived, the Landsat Science Team formed the basis for a continuing joint 
government and university working group for calibration and validation. 
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2. This Landsat calibration working group brought together the government IAS staff 
with other university based teams primarily involved in the vicarious calibration of 
satellite instruments. Regular redundant vicarious measurements of Landsat calibration 
had been missing up until this point. 

3. The continued NASA support for the Land cover Project Science Office (LPSO) 
allowed continued NASA participation in the post launch calibration and 
characterization mix, working alongside the IAS staff. LPSO funding covered the 
vicarious measurements beginning about 5 years after launch as well. 

D. Closer working relationship with national standards labs (NIST) 

Transfer radiometers directly calibrated at NIST were used to validate and adjust the 
radiometric calibration applied to the ETM+ sensor prior to launch (Markham et al,. 
1998). These efforts were part of the larger EOS calibration office efforts (Butler et al., 
2003). This provided better characterization of the linearity of the ETM+ as well as 
reduced uncertainty in the absolute radiometric scale. 

E. Pseudo Invariant Calibration Sites (PICS) 

The identification of pseudo-invariant sites by a number of investigators (e.g., Cosnefroy 
et al., 1996) provided an additional methodology for monitoring the stability of the 
Landsat-7 ETM+. The development of a long term acquisition plan (LTAP) for 
Landsat-7 significantly increased the number of times various sites were acquired 
around the world in the first year or so of Landsat-7 operations. After that point in time, 
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the identified PICS were given a higher priority in the LTAP leading to more frequent 
acquisition. 

2.2.2 Landsat-7 Performance over 10 years (Markham et al., 2011) 

The radiometric stability of the Landsat-7 ETM+ reflective bands has been continuously 
monitored over its 12+ year lifetime using the on-board calibration devices with the IAS as well 
as vicarious calibration and PICS (Markham et al, 2004; Markham et al, 2011). In general, the 
ETM+ appears to be more stable than any device or technique used to monitor it to date. The 
on-board lamp 1, which is on during the vast majority of imaging, shows a change in ETM+ 
response of nearly 30% in band 1 and decreasing changes in the longer wavelength bands 
(Table 3; Figure 1). Lamp 2, which is used about 0.1% of the time, shows changes about one 
quarter of those observed for lamp 1, so clearly lamp or the lamp optical system change is a 
significant contributor to the overall changes. The diffuser data shows overall changes similar 
in magnitude to the lamp 2 results, but spectrally different, with the largest decreases in the NIR 
region, very little change in the SWIR1 band and a small increase in the SWIR2 band. This is 
apparently the signature of the paint used on the solar diffuser. The vicarious measurements 
have the largest scatter of the various techniques and therefore are difficult to use to detect 
trends in response. The PICS sites, though also with significant scatter and some systematic 
annual-like variation, have the advantage of a larger sample of data points. Currently, the PICS 
sites appear to be the most useful for detecting long-term trends in ETM+. The apparent trends, 
some of which are statistically significant, though the systematic variations make the statistics 
suspect, are generally small, i.e. less than 1% or so in the VNIR bands and circa 2% in the 
SWIR bands total over 12 years. An independent study using PICS sites (Chander et al., 2010) 
showed changes in ETM+ response of 2% or less over the first 10 years of the Landsat-7 
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mission and trends generally consistent with MODIS. These numbers likely give reasonable 
bounds to the amount of time-dependent error being introduced by using a fixed calibration over 
the current mission lifetime. Work is ongoing to try to remove some of the systematic variation 
in these datasets, e.g. due to atmospheric and BRDF effects, to improve their utility for long- 
term instrument performance monitoring. 

The requirement on absolute calibration of the ETM+ was ±5% (one sigma). Comparisons of 
the operational calibration and the initial calibration of the ETM+ using the solar diffuser (based 
on the pre-launch measured reflectance of the diffuser and published solar irradiance curves) 
showed agreement to 4% or better overall and 2% or better in the VNIR (Markham et al, 2004). 
Estimates of the absolute calibration uncertainty using the vicarious measurements discussed 
above, i.e., the averages of the vicarious calibrations versus the operational calibration are 
presented in Table 4. Again, the differences are generally less than 3%, with the exception of 
the SDSU results in the two SWIR bands. Note, however, the vicarious measurements are 
almost all low relative to the operational calibration. Overall, results are consistent with the 
ETM+ radiometric calibration uncertainty being 5% or less, perhaps as good as 3% and can 
serve as a good standard for the Landsat series of sensors. 

2.3 Reconstruction of Landsat-5 TM calibration and tie to Landsat-7 ETM+: 2001- 
2007 

There are a number of aspects to this calibration reconstruction that include: 1) artifact 
correction, primarily memory effect and scan correlated shift to improve the internal 
consistency of the data (Helder et al, 2004), 2) characterization of the effects of the build-up of 
ice on the cold focal plane window which affects the calibration of TM bands 5 and 7 (Helder 
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and Micijevic, 2004) 3) development of the time history of the calibration based on on-board 
calibration sources, vicarious calibration and PICS (Teillet et al., 2004; Helder et al., 2008) and 
4) anchoring the time history of the calibration to Landsat-7 ETM+ (Teillet et al., 2001). A 
synopsis of the third and fourth aspects will be presented here. 

2.3.1 Phase 1: 1999- 2003 (Teillet et al., 2004) 

At the time of the Landsat-7 launch in 1999, a private company (EOS AT) had been operating 
Landsat-5 for approximately 13 years. The emphasis of the radiometric processing during this 
period of time and before was on providing visually clean products, i.e., free from striping and 
banding and other artifacts. With the launch of Landsat-7 and commencement of operations 
under US government control in 1999, reduced data prices and increased data acquisitions, there 
was renewed scientific interest in the Landsat data for long term studies. This was further 
enhanced with the return of Landsat-5 to US government control in 2002. This development led 
to a new, largely retrospective look at the Landsat-5 TM radiometric calibration. 

Initial work emphasized the cross calibration of Landsat-5 TM with the recently launched 
Landsat-7 ETM+. A cross calibration opportunity with Landsat-5 TM had been designed into 
the commissioning phase of Landsat-7. Lor three days, the Landsat-7 ETM+ imaged areas 
within about 15 minutes of imaging by Landsat-5 TM. The near simultaneous data sets, along 
with ground reflectance and atmospheric measurements made at several sites were used to 
perform the cross calibration (Teillet et al., 2001). 

In nearly the same time frame, reconstruction of the calibration record was begun. At this time, 
the information available was from: (1) pre-launch measurements, (2) the lamp-based internal 
calibrator system and (3) a relatively sparse record of vicarious calibrations performed by the 
University of Arizona. Given that the lamp and vicarious calibrations both showed some trends, 
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a fixed calibration, pre-launch or otherwise, was not appropriate, at least not for all bands. The 
internal calibration (IC) record, if limited to the best lamp, was smooth and discontinuity free. 
For this reason, a modified lamp record, using the best lamp and removing trends believed to be 
inherent in the calibration system itself was used as the basis for the calibration history (Teillet 
el al, 2004). 

The IC record for Landsat-5 when tied to an absolute scale via the June 1999 cross calibration 
with Landsat-7 ETM+, gave reasonably consistent results with the limited number of vicarious 
calibrations performed by the University of Arizona over the mission life (Figure 2), so this 
result was selected to represent the TM calibration lifetime. Termed “Look Up Table 2003” 
(LUT03), it was used by the operational system to process Landsat-5 TM data from May 2, 
2003 to April 20, 2007. 

2.3.2 Phase 2: circa 2007 (Helder et al., 2008) 

It was recognized when LUT03 was developed and released that it was a significant 
improvement over the previously used calibration technique, but that there was still some 
investigation needed. In particular, the vicarious calibration record was absent for a number of 
years (1988-1992), the vicarious calibration procedures had changed over the years and the 
selected LUT was not always a good fit to the existing vicarious data. In an attempt to fill this 
gap in vicarious data and to provide a more continuous and consistent calibration record, 
“pseudo-invariant calibration sites” (PICS) were investigated. 

Cosnefroy et al., (1996) identified a candidate list of desert sites useful for the monitoring the 
radiometric stability of satellite sensors. All these sites were in North Africa or on the Saudi 
Arabian peninsula. Unfortunately coverage of any of these sites with Landsat-5 TM data that 
was in the USGS archives was extremely limited, i.e., insufficient to develop trends. Coverage 
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of North America was quite robust, so a site in one of the driest areas of North America, the 
Sonoran desert in Mexico south of Yuma, Arizona was chosen. By treating some areas in the 
Sonoran desert as invariant sites, candidate responsivity trends were developed. Unfortunately 
these data had considerable scatter and the correct trends were not obvious. 

Searching the archive of Landsat-5 TM data revealed that path 181 row 40, essentially Libya 4 
in Cosnefroy’s terminology, was well represented in the European Landsat data archive. 
Through the cooperation of ESA, a good sampling of the Landsat-5 coverage of this site was 
obtained. These data, adjusted for illumination variation due to solar zenith angle and Earth- 
Sun Distance, showed clear trends, at least for TM bands 1-3 (Figure 2, for example). 
Independent work by de Vries et al (2007) showed similar results, so these trends, anchored to 
the ETM+ cross calibration point, formed the basis for LUT07. Effective April 21, 2007, 
LUT07 became operational in the US Landsat data processing system for Landsat-5 TM. 

2.4 Reconstruction of Landsat-4 TM calibration history and tie to Landsat-5 TM: 
circa 2008 (Helder et al, 2011a) 

Once the proper calibration for Landsat-5 TM had been developed, the next logical step was to 
concentrate on Landsat 4 TM. Even though Landsat 4 and 5 largely overlapped one another 
chronologically, calibration of Landsat-4 TM was desirable to extend the coverage from 1984 
back to 1982, and because Landsat 4 data complemented Landsat-5 data in that it was tasked to 
cover international land areas, whereas Landsat-5 focused more on the CONUS region. 

Initially, Landsat-4 TM was calibrated prior to launch through viewing of an integrating sphere 
that was traceable to the National Bureau of Standards (NBS, now National Institute for 
Standards and Technology or NIST), and this calibration was transferred to the onboard 
calibration system (which was identical to Landsat-5 TM). Following launch, the onboard 
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calibration system was used to calibrate all data collected by the instrument throughout its entire 
lifetime. An initial look at cross-calibration to Landsat-5 was conducted by Metzler and Malila 
(1985) during an under-flight opportunity when Landsat-5 was first placed into orbit. Their 
results indicated the two instrument gains were within 15% of each other. 

More recently, efforts were initiated to establish an accurate lifetime calibration record for 
Landsat-4 TM. These activities were begun by Mettler (2005) who used two PICS locations — 
the Sonoran Desert near Yuma, AZ, and a location near the Iraqi/Kuwait border to establish a 
cross-calibration with Landsat 5 through use of near simultaneous overpasses. These locations 
were driven primarily by the availability of data pairs within eight days of each other. Results 
from this work indicated that the reflective band calibration for Landsat-4 TM was within 5% of 
Landsat-5 TM in 1988 and 1991. An extension of this work by Malla (2008) using scene pairs 
that spanned dates from 1988 through 1993 produced similar results indicating the calibration 
for Landsat-4 TM had changed less than 5% over this extended time period. 

These results showed good stability for the TM sensor on Landsat 4, but were still limited to 
point measurements. Trends recorded by the onboard calibration system were very similar to 
what had been observed from Landsat-5 TM, did not agree with the results just mentioned, and 
were deemed not useful for a temporal calibration of the instrument. So, to establish a 
continuous calibration record for the instrument, Helder used PICS locations in Sonora and 
Libya 4 (Helder 201 la). An example of these results is shown in Figure 3. This plot shows the 
absolute calibration data (diamonds) that was obtained in 1988 and 1991, along with trending 
data from the two PICS shown as triangles. Even though the PICS trending data is sparse, there 
is a clear downward trend in the calibration of band 1 as shown by the black line. The 
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remaining reflective bands for Landsat-4 did not indicate any significant trend, and so the gain 
models for those bands are constant. 

2.5 Reconstruction of Landsat-1 to Landsat-5 MSS calibration history: circa 2010 
(Helder et al., 2011b) 

As mentioned previously, the first five Landsat satellites carried MSS sensors. These sensors 
acquired imagery of the earth’s surface from 1972 through 1992. As such, they add not only 
another decade of Landsat acquisitions, but also add to the depth and breadth of the historical 
Landsat archive. 

Calibration of MSS data was historically performed through use of an onboard calibration lamp 
system (Markham and Barker, 1987). Like TM, the MSS was a whiskbroom scanning 
instrument and, at the end of every second scan, the detectors were able to view the calibration 
lamp through a neutral density filter. Six different intensity levels from the lamp were recorded 
and regressed against coefficients derived during prelaunch calibration efforts to obtain gain and 
bias coefficients. Because the detectors in bands 1 through 3, which covered the visible and 
near infrared regions from 0.5 to 0.8 microns, were photomultiplier tubes, it was necessary to 
calculate multiple gain and bias coefficients during an image collect and combine them through 
use of a smoothing filter to minimize striping due to detector responsivity changes. An 
additional scene-based de-striping algorithm was implemented with Landsat 4 and 5 to equalize 
detector means and standard deviations. 

MSS data was stored in the USGS EROS archive in a variety of formats commonly known as 
MSS-X, MSS-P, and MSS-A (Helder, 1993). The MSS-X data format was used for Landsats 1- 
3. These data are stored in a radiometrically corrected format using the regression approach 
with onboard calibration lamps described previously. MSS-P data were archived for a short 
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period of time from 1979 through 1981. This data format stored the imagery corrected both 
radiometrically and geometrically. Lastly, the MSS-A format which was used for all of Landsat 
4 and 5 data was only radiometrically corrected. 

Since the historical calibration of MSS data was based on use of the onboard calibration lamp 
system with prelaunch gains and biases, with no updates during the lifetime of the instruments, 
significant differences in the calibration from one MSS sensor to another were known to exist. 
It was decided that cross-calibration of these sensors could best be established through the use 
of PICS as was done for Landsat-4 and -5 TM (Helder 2011b). The most difficult part of this 
activity was identification of appropriate PICS data sets. Since most of the MSS data available 
in the USGS EROS archive were of North America, the Sonora site was used extensively. 
Because of the potential for absolute calibration of the MSS sensor on Landsat-5, all of the other 
MSS sensors were cross-calibrated to MSS-5 in a stepwise fashion. That is, Landsat- 1 was 
cross-calibrated to Landsat-2; Landsat-2 was cross-calibrated to Landsat-3, and so forth through 
to the cross-calibration of Landsat-4 with Landsat-5. 

A second difficulty with the cross-calibration of the MSS sensor series was the significant 
differences that existed in the spectral bandpasses of the sensors. This is illustrated in Ligure 4 
which compares the average relative spectral response functions for MSS-3 to MSS-4. It was 
estimated that 1% error was possible simply due to spectral response differences. Lastly, since 
coincident collects were not available for most of the image pairs, the cross-calibration had to 
rely on near-coincident imaging events that occurred up to several days apart. This drawback 
allowed changes in the atmosphere to also affect the cross-calibration. Estimates of uncertainty 
are on the order of 4% for bands 1-3 and up to 12% for band 4 due to the presence of a large 
water absorption feature in this particular band. 
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Despite these difficulties, good cross-calibration results were achieved across all MSS sensors. 
An example is shown in Figure 5 for MSS-2 and MSS-3. In this figure the cross-calibration 
across all four spectral bands is shown using three different scene pairs from 1979 through 
1981. In each case, the imagery for each scene pair was obtained nine days apart. The data 
were very highly correlated and spanned a good portion of the dynamic range of both 
instruments. Unfortunately, a bias difference was discovered in three spectral bands for MSS-1 
through MSS-3. This can be seen for two of the three cases in Figure 5 for bands 2 and 3. In all 
other band/sensor combinations a simple gain factor was all that was necessary to place the 
instruments on a common radiometric scale. 

As an example of the differences that existed between MSS sensors, Figure 6 shows the band 1 
response of each sensor to the Sonora PICS. If the sensors were properly calibrated, each one 
would indicate the same radiance level for this invariant site. However, the figure clearly shows 
the average radiance level recorded by MSS-1 was greater than 154 radiance units, while the 
lowest level recorded was by MSS-4 at approximately 134 radiance units, a difference on the 
order of 15%. In addition, this figure also shows that the response of MSS-2 and MSS-3 
changed as a function of time as indicated by the dotted regression lines. Thus, a time 
dependent factor needed to be introduced to the calibration equation for these two bands (as 
well as for band 2 in MSS-2). 

After calibration has been applied to the data obtained from the Sonora PICS, the consistency of 
the sensors is substantially improved as shown in Figure 7. Now each sensor’s average 
response to the site is within 151.7 ± 0.05. Results for the other bands were similar with the 
exception of band 4, which still exhibits an uncertainty on the order of 5% due to atmospheric 
absorption features. 
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2.6 Tie of Landsat MSS to Landsat TM radiometric calibration: circa 2011 (Helder 


et al., 2011c) 

Absolute calibration of the MSS sensors was accomplished by taking advantage of the fact that 
Landsat 5 carried both an MSS and a TM sensor (Helder 2011c). This provided many 
opportunities for coincident collection of imagery. Since Landsat-5 TM was well calibrated, it 
was used as a transfer radiometer for the calibration of MSS-5. From there it was a simple 
matter to adjust the first four MSS sensors appropriately based on the results from the previous 
section. 

It should be emphasized that the TM and MSS sensors employ significantly differing spectral 
responses as clearly shown in Figure 8. Note especially that band 4 is much broader with the 
MSS sensor, and band 3 MSS really has no counterpart with respect to TM. Thus, the effort to 
place MSS on an absolute scale critically depends on an accurate knowledge of the RSRs of the 
two sensors, as well as the spectral signature of the test site. Secondly, it must be emphasized 
that MSS and TM cannot be expected to give the same radiometric response to targets with 
varying spectral signatures. To take into account the spectral differences between the two 
sensors, spectral band adjustment factors were calculated for each MSS band. These factors are 
basically the ratio of the radiometric response of the two sensors as shaped by their respected 
spectral filters in conjunction with the spectral signature of the target. 

Once again the Sonora Desert PICS location was used for the absolute calibration of MSS. This 
was due to the well-known stability of the site, invariance to slight amounts of misregistration, 
and also due to the availability of coincident data pairs. For this analysis, 25 scene pairs were 
utilized spanning a large portion of the lifetime of MSS 5. Several locations were used to span a 
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large dynamic range of the sensors. Lastly, three Hyperion scenes were available for the Sonora 
site and these served to guide the selection of the ROIs as well. 

Table 5 shows the gain ratios obtained for MSS-5 to TM-5. For all bands, the MSS response is 
slightly less than the TM response, ranging from 9% to 18%. Good repeatability was obtained 
with all estimates within 3% of each other. 


3. Reconstructed Calibration 

3.1 Using the data 

As a result of the efforts described in this paper, all of the Landsat sensors have been placed on 
a consistent absolute radiometric scale. Tables 6 through 10 provide the radiometric gain 
coefficients for all reflective bands of all Landsat sensors (excluding RBV) to date. Table 6 lists 
the gain coefficients for MSS 1 — 5 and can be applied according to 

La. = G * DN + B (1) 

with spectral radiance, La, given in units of W/(m -sr-pm) and gain and bias referenced by ‘G’ 
and ‘B’ respectively. In addition, for three bands (MSS-2 bands 1 and 2, and MSS-3, band 1) a 
time dependent factor (TDF) must be utilized as given in Table 7. For these bands, the TDF 
must be calculated as shown in the table as a function of time where T and TLaunch are in units of 
decimal years. The TDF is a multiplicative factor that is applied to both the gain and bias 
coefficients for those respective bands. Landsat 4 radiometric gains are given in Table 8. For 
this instrument, all gains are constant except for band 1. For band 1, radiometric gain is 
determined by multiplying the slope term (-0.0000418) by the number of days since launch and 
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adding the intercept term. Spectral radiance can then be calculated using the equation above. 
Sensor bias, B, is determined from onboard calibration shutter information at the end of each 
scan. Landsat 5 TM has an exponential model for all bands according to 

G„ e (p = a 0 e- a ' (t -‘°> +a 2 

where Gnew (t) is the time dependent radiometric gain, time, t, is in years and to = 1984.2082 
(Helder 2008). The gain model coefficients for Landsat 5 TM are listed in Table 9. Finally, 
Landsat 7 radiometric gains, Table 10, have not changed and can be used at any time. 

Fortunately for the Landsat data user, the above radiometric calibrations have all been applied to 
the Landsat data ordered from USGS EROS. Since Landsat data products are now available 
from USGS EROS free of charge, users are encouraged to reorder their data sets to obtain the 
most recent, and accurate, calibration. 

Data currently ordered from USGS EROS has had radiometric calibration applied to the 
standard terrain corrected product (LIT). The DN values associated with the pixels have been 
rescaled to a normalized dynamic range with a maximum radiance value of L M ax and a 
minimum value of Lmin (normally near zero). Converting the DN values in the LIT product 
involves a simple linear rescaling according to 
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Where Q ca imax and Q ca imin are 255 and 1, respectively, and LMAX>and LMIN^are part of the 
metadata that is included in the LIT product. 

3.2 Limitations of calibration and uncertainties 

Several of the papers cited previously (Helder 2008, Helder 2011a, b, c) form an extended 
calibration record for the Landsat series of sensors. Each of these papers gives the uncertainties 
for the cross-calibration of particular sensor pairs. However, no overall calibration uncertainty 
for the entire set of Landsat sensors has been presented. If the uncertainties for each sensor pair 
are considered to be uncorrelated with one another, then the total absolute uncertainty can be 
calculated through properly summing the individual components orthogonally. Table 11 
presents the results of these calculations; values have been rounded to integers, with no real loss 
of precision, for ease in observing trends. Here it is assumed that the uncertainties for the 
TM/ETM+ series of instruments are 5% for each instrument and each band. The uncertainty for 
the Landsat 5 TM to MSS calibration process is assumed to be similar to the MSS to MSS 
calibration process. Finally, for the MSS series, the uncertainties are propagated from sensor to 
sensor back to MSS-1. Pertinent observations from Table 11 indicate that all TM sensors 
exhibit progressively degraded calibration back to TM-4, which is within 9% absolute. 
Surprisingly, the MSS series of sensors is also within 10% absolute radiometric accuracy in all 
but band 4. Unfortunately, band 4 is very broad spectrally and includes large atmospheric 
absorption features. Thus, the uncertainty for this band in the MSS series grows much more 
rapidly to the point where the absolute uncertainty for MSS-1 band 4 approaches 25%. Landsat 
data users need to be aware of these limitations, especially for multi -temporal studies. 
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4. Summary 

Landsat sensors have recorded the Earth’s changes from 1972 through 2011 and provide the 
most complete remote sensing history of our planet available. Until recently, these sensors were 
not consistently calibrated and the ability to perform extended studies of the earth was limited. 
Fortunately, recent work has placed all Landsat radiometers (MSS, TM, and ETM+) onto a 
consistent radiometric scale, and radiometric uncertainties have been quantified. Coefficients 
for these radiometric gain models were presented. Better yet, these calibration updates have 
been instituted at the USGS EROS archive and all data currently being distributed contain the 
most accurate radiometric calibration available. Users of Landsat data now have an 
unparalleled opportunity to study the history of Earth’s land surfaces over the past 40 years. 
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7. Figure Captions 


Figure 1: Landsat-7 ETM+ Band 3 responses to on-board calibrators, vicarious calibrations and 
PIC sites over life of mission. 

Figure 2: Landsat-5 TM Band 2 responses to on-board calibrators, vicarious calibrations and 
PICS sites and the history of operational calibrations applied to the data. 

Figure 3. Absolute calibration model for Landsat 4 Thematic Mapper, band 1. 

Figure 4. Relative spectral response function of Landsat-3 and Landsat-4 MSS sensors. 

Figure 5. Landsat-2 to Landsat-3 cross-calibration for band 1 through band 4. 

Figure 6. Time series plot of TOA Radiance from Sonoran Desert site measured by band 2 of 
Landsat- 1 through -5 MSS sensors 

Ligure 7. Time series plots of TOA Radiance from the Sonoran Desert site as measured by band 
1 of Landsat- 1 through -5 MSS sensors after applying cross -calibration gains, biases and time 
dependent factor 
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8. Table Captions 


Table 1: Landsat Systems, Operational Dates and Sensors. 

Table 2: Landsat Sensor Characteristics. 

Table 3: ETM+ Apparent Changes in Gain based on various on-board calibrators, vicarious 
calibrations and pseudo-invariant sites. 

Table 4. Differences Between Operational ETM+ Calibration and Vicarious Calibrations: 
(Vicarious-Operational)/Operational. 

Table 5. Absolute Gain Ratios for MSS-5 to TM-5. 

Table 6. Radiometric Gain and Bias values for MSS 1 — 5. 

Table 7. Time Dependent Factors for MSS-2 and MSS-3. 

Table 8. Calibration Gain Model Coefficients for Landsat 4 TM. 

Table 9. Calibration Gain Model coefficients for Landsat 5 TM. 

Table 10. Gains for Landsat 7 ETM+. 

Table 11. Landsat Absolute Radiometric Calibration Uncertainties (%). 
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9.Figures 
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Figure 1. Landsat-7 ETM+ Band 3 responses to on-board calibrators, vicarious calibrations and PIC sites over life of mission. 
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Figure 2: Landsat-5 TM Band 2 responses to on-board calibrators, vicarious calibrations and PICS sites and the history of 


operational calibrations applied to the data. 
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Landsat 4 Band 1 Absolute Gain Model 



Figure 3. Absolute calibration model for Landsat 4 Thematic Mapper, band 1. 
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L3 and L4 MSS Band 1 RSR 



L3 and L4 MSS Band 2 RSR 



L3 and L4 MSS Band 3 RSR 


L3 and L4 MSS Band 4 RSR 




Figure 4. Relative spectral response functions of Landsat-3 and Landsat-4 MSS sensors. 
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Band-1 


y = 1.028x- 3.827 
R J = 0.997 


Band-2 


y=1.0700x -7.1890 
R- = 0.9967 



1x1 Line A Yuma-1 98 1(22-20) 

♦ Yuma-1 979(85-82) ♦ Yuma-1 979(38-35) L3/L2w/Btas 

L3/L2 w/o Bias 


Band-3 


y = 1.019x- 8.601 
R J =0.994 



1x1 Line A Yuma-1 98 1(22-20) 

♦ Yuma-1 979(85-82) ♦ Yuma-1 979(38-35) L3/L2w/B«as 

L3/L2 w/o Bias 


Band-4 


y=1.033x -1.179 
R : = 0.989 



1x1 Line A Yuma-1 98 1(22-20) 

♦ Yuma-1 979(85-82) ♦ Yuma-1 979(38-35) L3/L2w/Bias 

L3/L2 w/o Bias 



1x1 Line a Yuma-1 98 1(22-20) 

♦ Yuma-1 979(85-82) ♦ Yuma-1 979(38-35) L3/L2w/Bias 

L3/L2 w/o Bias 


Figure 5. Landsat-2 to Landsat-3 cross-calibration for band 1 through band 4. 
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Landsat 1-5 MSS Band 2 (TOA Radiance Vs Time) 
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Figure 6. Time series plot of TOA Radiance from Sonoran Desert site measured by band 2 of Landsat-1 through -5 MSS sensors 
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TOA Radiance derived over Sonoran Desert as apparently seen by Landsat-5 since 

1972, Band 1 
170 


> 160 

E 

| 150 

1 

7 140 

0 

c 

2 

1 130 

DC 

O 120 
110 
100 

1972 1977 1982 1987 1992 1997 

Digital Year 

♦ Landsat-1 ■ Landsat-2 * Landsat-3 »Landsat-4 x Landsat-5 







4 ♦ 




• 

♦ ■ 


X 

x * 

• 

X 
















L -1 

151.89 

L-Z 

151.69 

L -3 

151.64 

L -4 

151.58 

L -0 

151.74 


Figure 7. Time series plots of TOA Radiance from the Sonoran Desert site as measured by band 1 of Landsat-1 through -5 MSS 
sensors after applying cross-calibration gains, biases and time dependent factor 


38 


Relative Spectral Response Relative Spectral Response 


RSR Comparison RSR Comparison 




Figure 8. Relative Spectral Responses (RSRs) of Landsat 5 TM and Landsat 5 MSS. Also plotted is the spectral signature of the 
Sonora Desert PICS as obtained by Hyperion. 
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9. Tables 


Table 1: Landsat Systems, Operational Dates and Sensors 


Satellite 

Launch Date 

Decommissioning Date 

Sensors 

Landsat-1 

July 23, 1972 

January 6, 1978 

MSS (4 band) 

RBV (3 band) 

Landsat-2 

January 22, 1975 

February 5, 1982 

MSS (4 band) 

RBV (3 band) 

Landsat-3 

March 5, 1978 

March 31, 1983 

MSS (5 band) 

RBV (pan only) 

Landsat-4 

July 16, 1982 

June 2001 (lastTM data 

TM (7 band) 

MSS (4 band) 



transmitter failed 1993) 



Landsat-5 

March 1, 1984 

Operating 

TM (7 band) 

MSS (4 band) 

Landsat-6 

October 5, 1993 

Failed to achieve orbit 

ETM (8 band) 


Landsat-7 

April 15, 1999 

Operating 

ETM+ (8 band) 
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Table 2: Landsat Sensor Characteristics 


Sensor 

/Satellite 

Nominal Reflective 
Spectral Bands 

Spatial Resolution 
(Instantaneous Field 
of View) 

Normal Radiometric 
Quantization 

Radiometric Calibration 
Capabilities 

MSS/ 

500-600 nrn 

79 meters 

6-bit (compressed) 

Shutter with calibration 

Landsat 1-3 

600-700 nrn 

u 

u 

wedge neutral density filter 


700-800 nrn 

u 

u 

Solar calibrator 


800-1100 nrn 

u 

6-bit (linear) 


MSS/ 

500-600 nrn 

83 meters 

6-bit (compressed) 

Shutter with calibration 

Landsat 4-5 

600-700 nrn 

» 

u 

wedge neutral density filter; 


700-800 nrn 

a 

ii 



800-1100 nrn 

u 

6-bit (linear) 


TM/ 

450-520 nrn 

30 meters 

8-bit (linear) 

Shutter with lamp transfer 

Landsat 4-5 

520-600 nrn 



optics, 3 calibration lamps 


630-690 nrn 



with sequencer, radiance 


760-900 nrn 



feedback control 


1550-1750 nrn 





2080-2350 nrn 




ETM+/ 

450-515 nrn 

30 meters 

8-bit (linear) 

Shutter with lamp transfer 
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Landsat-7 


525-600 nm 
630-690 nm 
760-900 nm 
1550-1750 nm 
2080-2350 nm 
500-900 nm 


15 meters 


4 


12 


optics, 2 calibration lamps, 
constant voltage control, Full 
Aperture Solar Calibrator, 
Partial Aperture Solar 
Calibrator (ref] 



Table 3: ETM+ Apparent Changes in Gain (percent) based on various on-board calibrators, vicarious calibrations and pseudo- 


invariant sites. 


Band 

Lamp 1 

Lamp 2 

Diffuser 

PICS (linear fit) 





Delta 

Uncertainty 

1 

-28.3 

-6.8 

-7.1 

-1.1 

0.4 

2 

-24.8 

-4.6 

-7.0 

-0.3 

0.2 

3 

-19.5 

-4.0 

-8.5 

-0.2 

0.2 

4 

-3.8* 

+0.8* 

-11.8 

-1.3 

0.4 

5 

-7.2 

-0.8* 

-2.4 

-2.2 

0.3 

7 

-1.9* 

-0.9 

1.9 

-2.7 

0.5 

PAN 

-10.0* 

-2.5* 

-9.0 

-0.4 

0.4 


* not monotonic 
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Table 4: Differences Between Operational ETM+ Calibration and Vicarious Calibrations: (Vicarious-Operational)/Operational 


Band 

University of Arizona (n~95) 

South Dakota State University (n = 31) 


Mean Difference (%) 

Standard Deviation (%) 

Mean Difference (%) 

Standard Deviation (%) 

1 

-3.3% 

2.4% 

1.2% 

4.9% 

2 

-2.7% 

3.3% 

-0.2% 

6.1% 

3 

0.6% 

3.1% 

-0.6% 

7.0% 

4 

-0.6% 

2.8% 

-1.1% 

5.2% 

5 

-2.9% 

2.4% 

-7.2% 

5.6% 

7 

-2.2% 

2.9% 

-7.6% 

6.9% 

Pan 

- 

- 

- 

- 
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Table 5. Absolute Gain Ratios for MSS-5 to TM-5. 


MSS Band :TM 
Band 

1 : 2 

2 : 3 

3:4 

4:4 

ROI 1 

0.823 

0.914 

0.949 

0.951 

ROI 2 

0.830 

0.918 

0.950 

0.961 

ROI 3 

0.826 

0.909 

0.945 

0.951 

Average Gain 
Ratio 

0.824 

0.914 

0.948 

0.955 

Percent 

Uncertainty 

1.66% 

1.24% 

1.17% 

2.94% 
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Table 6. Radiometric Gain and Bias values for MSS 1 — 5. 


Landsat 1 MSS 

Landsat 2 MSS 

Landsat 3 MSS 


G 

2 

(W/m sr urn) 
DN 

Bias 

2 

(W/m sr |jm) 

G 

2 

(W/m sr urn) 
DN 

Bias 

2 

(W/m sr |jm) 

G 

2 

(W/m sr urn) 
DN 

Bias 

2 

(W/m sr \im) 

Band 1 

1.5968 

0 

*1.8036 

*7.1860 

*1.7507 

*3.4876 

Band 2 

1.2897 

9.1157 

*1.3150 

*0.7062 

1.2724 

2.7543 

Band 3 

1.3415 

-8.4567 

1.1520 

-2.4442 

1.0517 

2.9496 

Band 4 

1.2522 

0 

0.9654 

3.5493 

1.0349 

0.9505 

Landsat 4 MSS 

Landsat 5 MSS 



G 

2 

(W/m sr urn) 
DN 

Bias 

2 

(W/m sr |jm) 

G 

2 

(W/m sr urn) 
DN 

Bias 

2 

(W/m sr |jm) 

Band 1 

1 .7365 

3.7699 

1 .7345 

2.4937 

Band 2 

1 .2452 

3.9535 

1.2679 

2.7447 
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Band 3 

1 .0774 

4.9938 

1 .0693 

4.7483 

Band 4 

0.8717 

3.9538 

0.9025 

2.8653 
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Table 7. Time Dependent Factors for MSS-2 and MSS-3. 


Landsat 2 MSS 


Time Dependent Factor 

Band 1 

147.72 

(0.56709*(r-r iaM J c+ft 144.85) 

Band 2 

170.85 

(0.53916* (T-T Lau X+ h l6S.ll) 

Landsat 3 MSS 

Band 1 

151.55 

(1 .5 25 1* (7 — 7 L „ n \+ h 1 44. 1 0) 
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Table 8. Calibration Gain Model Coefficients for Landsat 4 TM. 


Model 

Coefficient 

Band 1 

Band 2 

Band 3 

Band4 

Band 5 

Band 7 

Slope 

(gain/DSL) 

-0.0000418 

0 

0 

0 

0 

0 

Intercept [gain] 

1.494 

0.719 

0.954 

1.073 

7.708 

14.65 
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Table 9. Calibration Gain Model coefficients for Landsat 5 TM. 



Model Parameters 

Band 

ao 

ai 

0C2 

1 

0.2901 

0.1399 

1.209 

2 

0.1246 

0.1045 

0.63 

3 

0.0839 

0.2386 

0.903 

4 

0 

0 

1.082 

5 

0 

0 

7.944 

7 

0 

0 

14.52 
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Table 10. Gains for Landsat 7 ETM+ 



Band 1 

Band 2 

Band 3 

Band 4 

Band 5 

Band 7 

Band 8 

High Gain 

1.225 

1.191 

1.538 

1.496 

7.589 

21.80 

1.483 

Low Gain 

0.8163 

0.7938 

1 .0245 

0.9969 

5.059 

14.532 

0.9885 
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Table 11. Landsat Absolute Radiometric Calibration Uncertainties (% 



L7 ETM+ 

L5 TM 

L4TM 

MSS 5 

MSS 4 

MSS 3 

MSS 2 

MSS 1 

Band 1 

5 

7 

9 

8 

9 

9 

10 

11 

Band 2 

5 

7 

9 

8 

9 

9 

10 

11 

Band 3 

5 

7 

9 

9 

10 

10 

11 

12 

Band 4 

5 

7 

9 

14 

18 

18 

22 

25 

Band 5 

5 

7 

9 






Band 7 

5 

7 

9 






Band 8 

5 
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